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ABSTRACT 

We use dark matter haloes identified in the MareNostrum Universe and galaxy 
groups identified in the Sloan Data Release 7 galaxy catalogue, to study the relation 
between halo shape and halo dynamics, parametrizing out the mass of the systems. A 
strong shape-dynamics, independent of mass, correlation is present in the simulation 
data, which we find it to be due to different halo formation times. Early formation time 
haloes are, at the present epoch, more spherical and have higher velocity dispersions 
than late forming-time haloes. The halo shape-dynamics correlation, albeit weaker, 
survives the projection in 2D (ie., among projected shape and 1-D velocity dispersion). 
A similar shape-dynamics correlation, independent of mass, is also found in the SDSS 
DR7 groups of galaxies and in order to investigate its cause we have tested and used, 
as a proxy of the group formation time, a concentration parameter. We have found, 
as in the case of the simulated haloes, that less concentrated groups, corresponding 
to late formation times, have lower velocity dispersions and higher elongations than 
groups with higher values of concentration, corresponding to early formation times. 

Key words: galaxies: clusters: general - galaxies: haloes - cosmology: dark matter - 
methods: N-body simulations. 



1 INTRODUCTION 

According to the hierarchical model of structure formation, 
groups and clusters of galaxies, embedded in dark matter 
haloes, emerge from Gaussian primordial density fluctua- 
tions and grow by accreting smaller structures, formed ear- 
lier, along anisotropic directions. Such structures constitute 
therefore an important step in the hierarchy of cosmic struc- 
ture formation and are extremely important in deciphering 
the processes of galaxy and structure formation. It should 
thus be expected that recently formed structures are highly 
elongated, reflecting exactly the anisotropic distribution of 
the large-scale structure from which they have accreted their 
mass, and dynamically young. Therefore the determination 
of the dynamical state of groups and clusters of galaxies and 
its evolution, which could be affected by a multitude of fac- 
tors - like environment, formation time, etc - is a fundamen- 
tal step in investigating the hierarchical galaxy formation 
theories. 



Both analytical calculations and N-body simulations 
have consistently shown that the virialization process will 
tend to sphericalize the initial anisotropic distribution of 
matter, and therefore the shape of haloes is an indication of 
their evolutionary stage. Furthermore, the halo shape, size, 
and velocity dispersion (among others) are important factors 
in determining halo member orbits and interaction rates, 
which are instrumental in understanding galaxy formation 
and evolution processes (eg., Jing & Suto 2002; Kasun & 
Evrard 2005; Avila-Reese et al. 2005; Bailin & Steinmetz 
2005; AUgood et al. 2006; Bett et al. 2007; Ragone-Figueroa 
& Plionis 2007; Maccio, Button, van den Bosch 2008 and ref- 
erences therein). It has been found that dark-matter haloes 
are triaxial but with a strongly preferred prolatness (eg., 
Plionis, Basilakos & Ragone-Figueroa 2006; Gottlober & 
Yepes 2007 and references therein) and that there is also a 
correlation between the orientation of their major axes and 
the surrounding structures. Such alignment effects, among 



relatively massive haloes, have been found to be particularly 
strong (eg., Faltenbacher et al. 2002; Kasun & Evrard 2005; 
Hopkins, Bachall & Bode 2005; Ragone-Figueroa & Plio- 
nis 2007; Faltenbacher et al. 2008 and references therein). 
There are indications that the mentioned correlation might 
arise, or strengthen, from re-arrangements of the orientation 
of the halo axes in the direction from which the last major 
merger event occurred (?). Hence, the shape of dark matter 
haloes might well be correlated with the large scale struc- 
ture, providing clues to the fashion in which mass is aggre- 
gated (eg., Basilakos et al. 2006). The imprints of this accre- 
tion can be observed in the substructure features present in a 
halo, which have been found to correlate with the halo shape 
(Ragone-Figueroa & Plionis 2007). Haloes with higher levels 
of substructure, or equivalently dynamically young haloes, 
are preferentially more elongated. The subsequent evolution 
should induce a change in the shape of haloes, leading to 
more spherical configurations, increasing their velocity dis- 
persions as well. 

According to this last statement, Faltenbacher & Math- 
ews (2007) have found, in agreement with the Jeans equa- 
tion, that the velocity dispersion of sub-haloes increases with 
the host halo concentration. Given that halo concentration 
and formation time are also correlated (Bullock et al. 2001; 
Wechsler et al. 2002), this result implies that older haloes 
should have higher sub-halo velocity dispersions. Once a 
halo is virialized, the dynamical mass can be computed from 
their velocity dispersion and virial radius (eg., Binney & 
Tremaine 1987). Nevertheless, if virialization has not yet 
been achieved, the above mentioned effect and its possible 
dependence to the local halo environment should be taken 
into account. In addition, there is also evidence for the spher- 
icalization of DM haloes with time (e.g. AUgood et al. 2006). 
This behavior can be well studied using haloes in numerical 
simulations, where their 3D shape and velocity dispersion 
can be calculated accurately. 

These issues motivated us to study the relation between 
the DM halo velocity dispersion and shape, in narrow halo 
mass ranges, in order to exclude the possible halo-mass de- 
pendent effects. We also investigate how this correlation is 
transformed when computing projected shapes and velocity 
dispersion along the line-of-sight, in order to be able to com- 
pare with observational data, like the SDSS galaxy groups. 

The plan of the paper is the following. In section [5] we 
describe the MareNostrum numerical simulation and the 
SDSS observational data together with the identification 
procedure of the corresponding systems and their proper- 
ties. In section [3] we study the halo shape-dynamics relation 
in the simulation. In section |4] we seek for a formation-age 
proxy which can be computed in a realistic group catalog. 
We continue in section [5] by presenting the shape-dynamics 
correlation results for the observed SDSS groups of galaxies. 
Finally, section |B] contains our conclusions. 



2 THE HALO & GROUP DATA 

2.1 The MareNostrum Simulation 

This non-radiative SPH simulation, dubbed The 
MareNostrum Universe (see Gottlober & Yepes 2007) 
and covering a volume of (500 h~^Mpc)^, was performed 



in 2005 with the parallel TREEPM-f SPH Gadget2 code 
(Springel 2005). The resolution of the simulation is such 
that the gas and the DM components are resolved by 
2 X 1024'^ particles. The initial conditions at redshift z = 40 
were calculated assuming a spatially flat concordance 
cosmological model with the following parameters: the total 
matter density Qm = 0.3, the baryon density Qb = 0.045, 
the cosmological constant Oa ~ 0.7, the Hubble parameter 
h = 0.7, the slope of the initial power spectrum n — 1 
and the normalization as = 0.9. This resulted in a mass of 
8.3 X lO^h'^MQ for the DM particles and 1.5 x lO^/i'^X© 
for the gas particles. The simulation followed the nonlinear 
evolution of structures in gas and dark matter from z = 40 
to the present epoch. Dissipative or radiative processes 
or star formation were not included. The spatial force 
resolution was set to an equivalent Plummer gravitational 
softening of 15 h~^ comoving kpc. The SPH smoothing 
length was set to the distance to the 40*'' nearest neighbor 
of each SPH particle. 

After the release of the 3-year WMAP data the original 
simulation has been repeated with lower resolution assum- 
ing the predicted low WMAP3 normalization (erg = 0.75) as 
well as with a higher normalization of erg = 0.8, which is bet- 
ter in agreement with the 5-year WMAP data. Yepes at al. 
(2007) argue that the low normalization cosmological model 
inferred from the 3 year WMAP data results is barely com- 
patible with the present epoch X-ray cluster abundances. All 
these simulations have been performed within the MareNos- 
trum Cosmology project at the Barcelona Supercomputer 
Center. 



2.1.1 Halo identification, shape & formation time 

In order to find all structures and substructures within the 
distribution of 2 billion particles and to determine their 
properties we have used a hierarchical friends-of-friends al- 
gorithm (Klypin et al. 1999, Gottlober et al. 2006b). At all 
redshifts we have used a basic linking length of 0.17 of the 
mean inter-particle separation to extract the FOF clusters, 
which corresponds to an overdensity of ~330 of the mean 
density. Shorter linking lengths have been used to study sub- 
structures. The FOF analysis has been performed indepen- 
dently over DM and gas particles. Using a linking length of 
0.17 at redshift z = we have identified more than 2 million 
objects with more than 20 DM particles which closely follow 
a Sheth-Tormen mass function (Gottlober et al. 2006a). 

Using the FOF method one extracts rather complex 
objects from the simulation which are characterized by an 
iso-density surface given by the linking length. In first ap- 
proximation these objects can be characterized by triaxial 
ellipsoids. The shape and orientation of the ellipsoids can 
be directly calculated as eigenvectors of the inertia tensor of 
the given object. Then the shape is characterized by the ra- 
tios between the lengths of the three main axes {a > b > c). 
Gottlober & Yepes 2007 have studied the shape of the galaxy 
clusters in the The MareNostrum Universe and found 
that both the gas and dark matter components tend to 
be prolate although the gas is much more spherically dis- 
tributed. Here we study the shape of the more massive dark 
matter component of the FOF objects which refiects the 
shape of the potential in which the galaxies move. 

Finally, in our present work we will extensively use the 
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Figure 1. Different panels correspond to tiie formation time 
distribution of haloes in a different range of mass, namely 1.0—1.3, 
2.0-2.5, 3.2 -4.0 and 6.3- 7.9 X lO^^h-^Ms, from left to right 
and top to bottom. Vertical dashed lines denote the first and 
fourth quartile of the distribution. 



notion of halo formation time (zform), which is usually de- 
fined as the redshift at which the halo accretes half of its 
final {z = 0) mass. By selecting haloes in the lower and up- 
per quartiles of the formation time distribution, we define 
what we will call the late formation time (LFT) and early 
formation time (EFT) halo populations, respectively. We 
also divide the haloes in four mass ranges, namely 1.0 — 1.3, 
2.0 - 2.5, 3.2 - 4.0 and 6.3 - 7.9 x lO^^h'^Me, since we 
wish to cancel the dependence of any halo-based parameter 
on its mass. 

Figure [T] shows the distribution of formation times for 
the four subsamples of haloes as defined previously. The ver- 
tical dashed lines denote the location of the first and fourth 
quartile of the distribution. 



2.2 SDSS DR7 Galaxy Groups 

In recent years a considerable effort has been invested in 
identifying, using a multitude of methods, groups and clus- 
ters of galaxies in redshift surveys (e.g. Huchra & Geller 
1982; TuUy 1987; Nolthenius & White 1987; Ramella et al. 
2002; Merchan & Zandivarez 2002, 2005; Gal et al. 2003; 
Gerke et al. 2004; Lee et al. 2004; Lopes et al. 2004; Eke 
et al. 2004; Tago et al. 2006, 2008; Berlind et al. 2006; 
Crook et al. 2007; Yang et al. 2008; Sharma & Johnston 
2009; Finoguenov et al. 2009). 

In this work we identify groups of galaxies in the SDSS 
DR7 spectroscopic galaxy sample, which comprises of more 
than 900000 galaxies, within an area of approximately 9000 
square degrees, and with a limiting r-band magnitude of 
JTiiim ~ 17.77. The group finding algorithm is based on 
the procedure detailed in Merchan & Zandivarez (2005) and 
consists in using the friends-of-friends algorithm, similar to 
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Figure 2. Fractional distribution of virial radii (left panel) and 
masses (right panel) of groups in a volume limited sample (with 10 
or more members) extracted from the SDSS DR7 corresponding 
to 0.06 < z < ziijjj = 0.10 (dashed histogram), compared to 
that corresponding to the total sample of groups (with 3 or more 
members) up to z = 0.10 (empty histogram). 



that developed by Huchra & Geller (1982). The algorithm 
links galaxies (i and j) which satisfy Dij < DoR{z) and 
Vij < VoR{z), where Dij is their projected distance and Vij 
is their line-of-sight velocity difference. The scaling factor 
R{z) is introduced in order to take into account the decre- 
ment of the galaxy number density due to the apparent mag- 
nitude limit cutoff. We have adopted a transverse linking 
length Do corresponding to an overdensity of 5p/p = 80 and 
a line-of-sight linking length of Vb = 200km s"'^. 

In order to avoid significant discreetness effects in the 
group shape and dynamics determination, we limit our 
catalogue to those groups with more than ten members 
(rim > 10). To further avoid artificial redshift-dependent 
effects (eg., Frederic 1995; Plionis et al. 2004; Plionis et al. 
2006) and to build a roughly volume limited sample, we se- 
lect those groups with 0.06 < z < 0.1 and Alio < Mun^, 
where Mio is the absolute magnitude of the lO"^ brighter 
galaxy and Miim(= —19.53) is the absolute magnitude of a 
galaxy with apparent magnitude m = mum at z = 0.1. The 
resulting sample, within the previously mentioned z-range, 
comprises 760 groups with at least ten member galaxies. 

For these groups, we compute their virial radius, line of 
sight velocity dispersion (ctz) and virial mass {A4), following 
Merchan & Zandivarez (2005). In Figure [2] we present the 
distribution of virial masses and radii of the resulting groups, 
comparing with those of the whole {rim > 3 and z < 0.1) 
parent group sample. It is evident that we are preferentially 
selecting the richest systems, which is the price we have to 
pay for producing a volume limited subsample. 



2.2.1 Group Projected Shape 

The group elongation is determined by diagonalizing the 
"inertia" tensor, which we construct weighting each galaxy n 
by its luminosity L„ . This is done in order to weight more the 
galaxies that dominate in shaping the group gravitational 
potential: 
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Figure 3. Left panel: The q — Qran scatter plot for groups within 
the volume limited sample and with at least 20 member galax- 
ies. Right panel: The median value of q — qran as a function of 
minimum number of group galaxy members. Errorbar correspond 
to the 33% and 67% quantiles. The solid line corresponds to the 
value q — qran = 0.063 and the dashed lines to the ±0.035 range. 



where Xi^„ is the i component of the position vector of 
the n"* galaxy relative to the group center of mass. Diago- 
nalizing the lij tensor we obtain two eigenvalues which are 
related to the square root of the minor and major axes, b 
and a respectively, which define the halo axial ratio: q — b/a. 

Furthermore, as it is well documented, there exists a 
dependence of the group halo apparent elongation on the 
member- number resolution of the system (eg., Paz et al., 
2006; Plionis, Basilakos & Ragone-Figueroa, 2006). Accord- 
ing to this effect, a given system will appear artificially more 
elongated (smaller q = b/a ratio) when sampled by a lower 
number of members. This effect however is large for groups 
with less than ten members (see fig. 3 in Plionis et al. 2006) 
and therefore it should not affect significantly our group 
sample, which by construction has rim > 10. In any case 
and to overcome any residual effect which could be intro- 
duced due to the variable resolution of the different groups, 
we sample each group having n™ > 10 by selecting randomly 
only ten galaxies, which is the number resolution of the low- 
est Um galaxy groups in our sample. Finally, the adopted 
value of the projected shape, qran, is the result of com- 
puting the mean of 10 realizations of the random galaxy 
selection procedure. By this procedure we degrade the accu- 
racy of the shape parameter of groups with rim > 10, but we 
gain resolution consistency over all our group sample. To get 
an idea of the uncertainty introduced by our procedure, we 
present in the left panel of Figure[3]the g vs. qran correlation. 
As expected, the qran values are systematically lower than 
their counterparts, q, computed using all member galaxies, 
while the deviation increases for the intrinsically more spher- 
ical systems. Since in this plot we show together all groups 
with Um > 20, in the right panel of FigOwe plot the me- 
dian value of 5 — qran as a function of minimum number 
of group member galaxies. It can be seen than there is, as 
expected, a tendency of g — qran to slowly increase with n,n 
but, within the errors, we can quote an overall median value 
of: q — qran — 0.063 ± 0.035, independent of n„i. 



3 THE HALO SHAPE-DYNAMICS RELATION 

In this section we investigate the relation between morphol- 
ogy and dynamical status of DM haloes. 

As it is well known halo properties can have a strong de- 
pendence on halo mass. Especially, the velocity dispersion 
of virialized or nearly virialized haloes are related to the 
halo mass through the virial theorem: M <x a'^ . Therefore, 
in order to search for an intrinsic dependence of the veloc- 
ity dispersion on shape, it is imperative to suppress their 
possible dependence to mass, due to the virialization expec- 
tation. To this end we have studied the shape-cr„ correlation, 
in the four mass ranges defined in section (2.1.2), both in 
3D-space and in projected 2D-space (using ID velocities and 
projected halo shapes). 



3.1 3D Halo c/a — o-u Correlation 

Figure |4] shows the isodensity contours of the shape-Uu scat- 
ter plot, at redshift 2 = 0, for some of the selected mass bins 
(left boxes of each panel). As it can be seen, there exists a 
clear correlation between halo sphericity and halo velocity 
dispersion, according to which more spherical haloes have 
higher velocity dispersions. Even though we have been cau- 
tious in selecting small mass ranges, the observed behavior 
could, in principle, be attributed to the halo mass variation 
within each mass-bin. For this to be true the more massive 
systems within each bin - which should also be the higher ve- 
locity dispersion ones - should have the highest values of the 
c/a ratio. However, it is well established that more massive 
systems tend to be less spherical (e.g. AUgood et al. 2006). 
Therefore, the observed correlation cannot be attributed to 
a possible residual dependence on halo-mass. A rather in- 
teresting possibility is related to the fact that haloes, of any 
given mass and at any given cosmic epoch, could have a 
different formation time as well as a different evolutionary 
history. 

Under this hypothesis, and in order to further investi- 
gate the dependence of shape-cr„ correlation on the halo for- 
mation time (zform), we plot separately the EFT and LET 
c/a — o-u correlations in the right boxes of each panel of 
Figure 3] It is evident that two populations occupy clearly 
delineated regions in the c/a — cr„ plane, with the EFT haloes 
(solid line) typically having higher velocity dispersions and 
higher sphericities than the LET haloes (dotted line) . It can 
be also observed that the LET contours cover a wider range 
of cr„ values than the corresponding EFT ones, but they 
do not reach the highest EFT sphericities. Finally, it is in- 
teresting to also note that there is a shift of the peak of 
the LET contours towards lower c/a's and (j„'s between the 
lower mass (upper panels) and the higher mass LET halos 
(lower panels) . This could be attributed to a fraction of the 
most massive LET systems being relatively more dynami- 
cally young and in the phase of primary mass aggregation. 
We will return to this issue later on. 

As it is obvious from Figure |4j the same overall trend, 
between the LET and EFT halos, is present in all mass bins. 
It would be ideal then, in order to enhance the statistical va- 
lidity of our results, to combine the information contained 
in the different mass ranges in a unique, independent of halo 
mass, analysis. This would also facilitate the corresponding 
study of the observational data, since the available number 
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Figure 4. Each panel corresponds to haloes in the four different ranges of mass, increasing in mass for the top left, to top right, then 
bottom left and bottom right panels, respectively. The left box of each panel shows the isodensity contours of c/a vs (t„ for all haloes in 
the corresponding range of mass. The right box of each panel shows the corresponding isodensity contours but only of the 25% and 75% 
quantiles of the formation time distribution, which correspond to the EFT (solid line) and LFT (dotted line) families of haloes. 



of SDSS groups is significantly smaller than that of the sim- 
ulated haloes and thus we cannot afford to divide them in 
many mass subsamples. 

Since the parameters under study (ct^, c/a and Zform) 
depend on halo mass, we have devised a normahzation proce- 
dure that imposes the mass independence of the results and 
combines in a unique relation the halo velocity dispersion 
and shape of the whole sample. In detail we derive normal- 
ization functions of the different parameters by computing 
analytic fits of the relation between the corresponding pa- 
rameter (P) and the halo mass {M). Then for a given halo 
mass we normalize each of the three parameter under study 
to the value expected from the analytic fit {P{M)). 

Since it is no longer necessary to split haloes in nar- 
row mass bins we can use the entire main sample of haloes 
to plot in Figure [5] the normalized, independent of mass, 
shape-dynamics correlation. It is evident that the EFT and 
LFT behavior, seen in Figure |3J is reproduced for the whole 
halo sample and we have further verified that these results 
are indeed independent of halo mass. In the inset of Fig- 
ure [5] we show that the corresponding distribution of masses 
of the EFT (solid line) and LFT (dashed line) populations 
are statistically equivalent and therefore there is no residual 



Table 1. The Spearman coefficient (SC) for the velocity 
dispersion-shape correlation for haloes in four mass bins. In all 
cases (3D and 2D) the correlations of the EFT haloes is signifi- 
cantly lower than that of the LFT haloes. 



Sample 



1.0-1.3 



Mass [xlO"X0/i-i] 
2.0-2.5 3.2-4.0 



6.3-7.9 



3D All 


0.37 


0.42 


0.44 


0.48 


3D LFT 


0.45 


0.55 


0.60 


0.72 


3D EFT 


0.11 


0.16 


0.21 


0.22 


2D All 


0.27 


0.31 


0.33 


0.34 


2D LFT 


0.31 


0.38 


0.42 


0.47 


2D EFT 


0.11 


0.12 


0.15 


0.18 



mass-dependent effect affecting the correlation shown in the 
main panel. 

Our main results, as extracted from Figures 4 and 5, 
can be tabulate as follows: (a) The LFT haloes show a sig- 
nificantly stronger c/a — a^ correlation than the EFT haloes, 
while they also span a wider range of values, (b) The LFT 
cr„ values are, in some cases, even higher than those of the 
EFT ("virialized") halo population, (c) A bimodal pattern 
is present in the (j„ — c/a plane (Figure [S} with two local 
maxima, one of which is situated at the low velocity disper- 
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Figure 5. Normalized c/a-cry relation for EFT (solid line) and 
LFT (dashed line) populations. All masses are considered. The 
formation time used to divide EFT and LFT is also normalized 
to the expected value for a given mass (see text). Inset panel: 
Distribution of masses of the EFT (solid line) and LFT (dashed 
line) populations demonstrating that no mass effect is present in 
the result obtained in the main panel. 



sion and low c/a values, which as we have discussed earlier 
and deduced from Figure 14] it is caused by the more massive 
haloes. Below we attempt to interpret individually each of 
the above results: 

(a) Different a^ — c/a correlation strength of the EFT 
and LFT haloes: The apparent lack of a strong correlation 
of the EFT halo population should be attributed to the fact 
that EFT haloes had plenty of time to virialize and that 
the virialization process erases the intrinsic shape-(j„ corre- 
lation. To quantify the strength of the shape-cr,j correlation 
we estimate the Spearman correlation coefficient as a func- 
tion of halo mass (Table [l]). Clearly the observed strong and 
significant correlation of the whole halo population should 
be attributed to the LFT haloes, since (as we already dis- 
cussed) the EFT haloes show very weak shape-(T„ correla- 
tions. 

So far, the results obtained suggest that the intrinsic 
halo shape-a„ correlation, within each small mass range, is 
the consequence of the coexistence of haloes of different for- 
mation times, going through different evolutionary phase^. 
Dark matter haloes, in the initial stages of formation and 
before virialization is complete, have typically a smaller ve- 
locity dispersion and a more elongated shape. Once virial 
equilibrium is reached these quantities stabilize and, since 
this equilibrium configuration is more likely to be achieved 
by the early formation time haloes, the EFT family shows 
only a weak such shape-(T„ correlation. 



Note that a similar conclusion was reached in a recent observa- 
tional study of 2^ 0.04 groups of galaxies (Tovmassian & Plionis 
2009). 
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Figure 6. Normalized c/a — cry trend corresponding to the 
entire sample of haloes. Each point corresponds to the median 
c/a/{{c/a){M)) and ay/{a\j{M)) values computed for haloes 
with 2form within different octiles of formation time. Error bars 
mark the 10% percentile of the corresponding distribution. 



We attempt to investigate in more detail the evolution 
paradigm, as the cause of the a"„ — c/a correlation, by di- 
viding the halo formation-time distribution in eight quan- 
tifies. We then compute the median values of a-y / {a{M)) and 
c/a/{{c/a){M)) within each Zform quantile and plot them in 
Figure [S] In concordance with Figure O we consistently find 
a clear trend indicating that the earlier the halo is formed, 
the higher its (t„ and c/a ratio, as expected from the fact 
that older haloes have more time to evolve and virialize. Ev- 
idently, a halo evolutionary sequence appears, with dynam- 
ically young haloes situated at the lower-left of the c/a — Oy 
plane and as they evolve, they move towards higher (t„ and 
c/a values, where the virialized haloes are located. 

Note, however, that in Figure [6] the dynamically 
youngest haloes (1st octile of formation time) show a be- 
haviour that deviates from the overall trend, being charac- 
terized by a lower normalized c/a value. This outlier has a 
relatively larger uncertainty with respect to the other oc- 
tiles, and a possible cause of this behaviour is presented in 
(c). 

(b) LFT haloes with o"„ values even larger than those 
of EFT haloes: In an attempt to explain what is the cause 
of the rather unexpected behavior of some LFT haloes, i.e., 
those having relatively high sphericity and in the same time 
a velocity dispersion which is as large or even larger than 
what expected from the "virialized" EFT haloes, we suggest 
two mechanisms: 

• Since LFT haloes are prone to have had a recent merger 
and since it is expected that dynamical interactions and 
merging processes increase the velocity dispersion of the in- 
volved systems some merging or highly interacting haloes 
could be undergoing a transient geometrical configuration of 
high sphericity (resulting when the in-falling substructures 
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Figure 7. Same as Figure \5\ but only considering haloes with 
masses between 3.2 and 7.9 X 10^^h~^A4Q. In the left panel the 
entire sample is used whereas in the right panel all double-core 
haloes have been excluded. 



are near the perigee), while temporarily achieving a velocity 
dispersion which could be even larger than their virial value 
(Faltenbacher et al. 2006). We have verified that indeed this 
is the case for some of the LFT haloes. 

• Part of the high velocity dispersion and quasi-spherical 
LFT haloes could be due to truly virialized LFT haloes by 
a somewhat rapid virialization process depending on their 
particular formation history and/or environment. 

(c) The two maxima in the LFT (t„ — c/a correlation 
plane: As we have discussed earlier, the lower a„ — c/a max- 
imum seen in Figure 5 is caused by the most massive haloes 
(ie., those with M^ 3.2 x 10^^ h~^ Mq) and could be at- 
tributed to a fraction of the most massive LFT systems being 
relatively dynamically younger and in the phase of merging, 
where the FOF algorithm could join a merging pair into a 
single FOF halo. 

In order to verify our suspision we have selected haloes 
with masses between 3.2 — 7.9h~ Mq (including the two 
most massive ranges presented in Figure U where the bi- 
modal pattern was evident for the LFT haloes) and inves- 
tigated in which manner the described trend would be af- 
fected if double (or multiple)-core halos (or otherwise dumb- 
bell shaped haloes) were excluded. Double-core haloes were 
identified as those haloes that are split in two or more sub- 
haloes of at least 1000 particles, within a relative distance 
of 2rsph (rsph is the diameter of a sphere which has the same 
volume as the considered halo) , when a re-identification with 
a shorter linking length (0.135 of the mean interparicle sep- 
aration) is performed. Therefore, they are major pre-merger 
systems which the FOF joins into a single system and which 
other halo-finders would have probably singled out each 
component as an individual halo. In the left panel of Figure 
[7] we show the c/a — a^ correlation when all haloes in the 
mentioned mass range are considered, whereas in the right 
panel double-core haloes are excluded. It is evident that the 
secondary maximum at low a^ and c/a, is caused by the 
double-core haloes, which are those that are in an active 
major merger process, and thus an integral part of the "dy- 
namically young" family of haloes. Note, that although the 
c/a — Ov LFT correlation is weakend when excluding the 
double (multiple)-core haloes, it is still clearly present. 

Returning to the issue of the "outlier" seen in Figure [H 
we have verified that the cause of the significantly smaller 




EFT Haloes 
LFT Haloes 



q / qW 

Figure 8. As in Fig[5]but using projected halo shapes and ID 
velocity dispersions. 



normalized {c/a) value of the first formation time octile 
is the significantly larger number of double (multiple)-core 
haloes found in this octile with respect to the other octiles. 
For example, for the case of the most massive haloes (ie., 
those investigated in Figure O we find that ~ 22% of the 
haloes in the first octile are double (multiple)-core haloes, 
while this number drops to 4.8% and 2.6% in the second and 
third octile respectively, and to even smaller values there- 
after. 



3.2 2D Halo q- a. Correlation 

Here we present the corresponding analysis of the previous 
subsection, but using the projected halo axial-ratio [q — b/a 
obtained from the projected particle distribution) and the 
one-dimensional velocity dispersion. This is done in order to 
see how the previous (3D) results are translated into theo- 
retical expectations that can be compared with the results 
based on observational group samples. 

Figure [8] is the 2D analogue of Figure [5] and we can 
clearly see that qualitatively the correlation revealed in Fig- 
ureO albeit weaker, is repeated in Figure[S] Details and sig- 
nificances may differ (see Table [l] for different mass ranges), 
however, the main result is that we expect to see (in realistic 
group samples) hints of an independent of mass correlation 
of the projected halo axial-ratio with its line-of-sight velocity 
dispersion. 



4 A FORMATION TIME PROXY 

The aim of this section is to provide a formation-age proxy, 
tested with simulated haloes, but that it can also be com- 
puted in realistic groups of galaxy samples, where the 
adopted definition of Zform is not usable. This is done in 
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Figure 9. The A^ubh parameter as a function of formation time 
computed for dark matter lialoes in four ranges of halo mass. 



order to test the morphology-dynamics-formation time cor- 
relation (which was found to be present in simulated haloes) 
in groups extracted from the SDSS DR7. 

We propose as a halo formation-age proxy (Asubh), the 
fraction of mass that is associated to the most massive sub- 
structure (subhalo), when a further re-identification is per- 
formed with a linking length which equals half the original 
one (0.17/2 of the mean inter-particle separation). This pro- 
posed parameter can be interpreted as a concentration in- 
dicator since it accounts for the distribution of mass at two 
different overdensities (given by the used linking lengths). 

The number of particles of a halo in each overdensity 
level is defined as rihigh and niow for the higher and lower 
overdensities, respectively. Then the concentration parame- 
ter is given by: 



Asubh = 



Mhigh 



(2) 



The parameter Agubh takes values within the range Agubh G 
[0, 1], with a value near 1 indicating a high concentration and 
a value near indicating the opposite. We show in Figure [5] 
the relation between Asubh and halo formation redshift for 
the different ranges of mass labeled in the plot. It is evident 
that the oldest haloes have the highest levels of concentra- 
tion, whereas the youngest ones extend to lower values of 

Asubh- 

Since the above analysis clearly demonstrates that the 
proposed parameter can be used as an age indicator, we 
continue with the computation of the concentration param- 
eter corresponding to galaxy groups in the observational 
data, Asub- To this end we perform a second friends-of- 
friends identification, but this time using an overdensity of 
Sp/p — 300. This procedure allows us to extract a denser 
subgroup, of a given system, with respect to the original 
Sp/p = 80 sample. Then Asub is given by the ratio between 
the number of members of the denser subgroup to the num- 



ber of members of its parent {Sp/p = 80) group. Utilizing 
this formation-age proxy, we examine the SDSS-DR7 group 
shape-dynamics-formation time relation in the next section. 



5 THE SHAPE- VELOCITY DISPERSION 
RELATION OF SDSS DR7 GROUPS 

Many studies have attempted to determine the morphol- 
ogy and dynamical status of groups of galaxies (eg. Hick- 
son et al. 1984; Malykh & Orlov 1986; Orlov, Petrova, 
Tarantaev 2001; Kelm & Focardi 2004, Plionis, Basilakos 

6 Tovmassian 2004; Plionis, Basilakos, & Ragone-Figueroa 
2006; Da Rocha, Ziegler, & Mendes de Oliveira 2007; Co- 
ziol & Plauchi-Frayn 2007; Wang et al. 2008; Hou et al. 
2009) and have found that their morphology corresponds 
to that of mostly prolate-like triaxial ellipsoids (see how- 
ever Robotham, Phillips & de Propris 2007), while there 
also appears to be a correlation between velocity dispersion 
and projected group shape (eg. Tovmassian & Plionis 2009) 
similar to what we found in the previous section. 

Here we analyse our "volume-limited" group catalogue 
based on the SDSS DR7 galaxy survey, which was described 
in section (2] 

As we have concluded in section [3^ the shape-cr corre- 
lation for dark matter haloes, extracted from the simulation, 
survives in projection but it is significantly weakened when 
the projected shape (q) and the one dimensional velocity 
dispersion are used. Further scatter should be expected in 
a more realistic situation in which the projected shape of 
SDSS groups are computed using only ten members (gran)- 
Indeed, as it was shown in section 12.21 the typical error, 
expected when only ten members are considered, in the es- 
timation of the projected shape is non-negligible. Such issues 
may act as to hide the morphology-dynamics correlation of 
observed groups of galaxies, and thus even a weak such cor- 
relation should be considered as a hint of a true underlying 
effect. We show in Figure [TO] the least square best fit line of 
the Qran VS. (Jz/(Tz{M) Scatter plot (of the entire sample of 
galaxy groups) where the observed trend is in full agreement 
with Figure [S] 

We now attempt to investigate whether the observed 
Qran VS. o , Correlation, in Figure 1101 is due to the coex- 
istence of groups of different formation times, as was in- 
deed shown to be the case in the dark matter halo analy- 
sis. As shown in section 3.2 (see also Figure [8} , using pro- 
jected shapes and one dimensional velocity dispersions does 
not completely mask the systematic behaviour observed for 
haloes of different formation times. However, taking into ac- 
count the unfeasibility of determining the formation time of 
the observed galaxy groups, we use as a proxy of the group 
dynamical age its concentration parameter, Asub- This pa- 
rameter has indeed been shown (for the case of DM haloes) 
to correlate with halo formation time (see section 3}- 

In Figure [10] we separately plot groups of the first and 
last quartile of the Asub parameter distribution (as was done 
for formation time in section [3TT} . Filled and empty circles 
represent galaxy groups with low and high concentration 
(corresponding to groups in the first and last quartiles of 
the Asub distribution), respectively. The tendency for high 
concentrated groups (open circles) to have systematically 
higher velocity dispersions and q values is quite evident. We 
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Figure 10. Qran vs. normalized a^ relation for groups in the 
SDSS DR7. Filled (low concentration) and empty (high concen- 
tration) dots stand for groups in the first and last quartile of the 
normalized Agub distribution, respectively. The large circles cor- 
respond to the median values. Solid line shows the least square fit 
using all the groups, irrespective of their Ag^b value. The inside 
box shows the mass distributions for groups in the cuartiles of 
low and high concentration (dotted and solid lines respectivelly) . 



also plot their corresponding median values as large filled 
and empty circles, from which it can be seen that indeed 
there is a separation of the two families (based on Asub), 
in both qran and a^ axes. We have again verified that the 
observed behaviour is not due to a mass-dependent effect, 
since the distribution of inass of the two group families are 
statistically equivalent (see insert of Figure [10]) . 

It is important to keep in mind that the uncertainty in 
the determination of the projected group shape, induced by 
our random sampling procedure aiming in reducing the vari- 
able resolution effects, could introduce a large scatter in the 
morphology-dynamics correlation. This could then act in the 
direction of hiding the differences between the two families 
of groups, represented by the different symbols in Fia llOl We 
have tested this assertion by using, within small ranges of 
group mass and for a given number of galaxy members, all 
the available galaxy members to estimate more accurately 
the group q values, in the expense however of working with 
a low number of groups. We now find significantly larger 
differences between the mean g-values of the low and high 
concentration families. Therefore, had we higher resolution 
measurements of group shapes, we would have recovered a 
significantly larger separation between the low and high con- 
centrated group families in the q — a^ plane, strengthening 
the interpretation of group formation time differences being 
the cause of the observed q — a, correlation, as is the case 
in the dark matter haloes analysed. 



6 SUMMARY & CONCLUSIONS 

We have identified dark matter haloes in the MareNos- 
TRUM Universe and groups of galaxies in the SDSS DR7 
with the purpose of searching for possible, independnet of 
mass, dependencies of halo/galaxy groups shape on their 
dynamics. 

We have computed the halo formation time as the red- 
shift at which the halo accretes half of its final mass. Since 
this definition of formation time is not applicable to observed 
groups identified in realistic galaxy catalogues, we have iden- 
tified a concentration parameter that can be used as a proxy 
of the formation time of haloes and at the same time it is 
easily computed for the observed SDSS galaxy systems. 

On this regard we have found a significant correla- 
tion between halo shape and dynamics (velocity dispersion) 
which is independent of the halo mass. We have identi- 
fied the cause of this correlation to be the halo formation 
time differences. Early formation-time haloes show a very 
weak such correlation, being more spherical, having a higher 
velocity dispersion and showing substantial concentration. 
On the contrary, late formation-time haloes show a strong 
shape-dynamics correlation (having a wide range of shape 
and velocity dispersion values) , with typically lower velocity 
dispersions, lower sphericities and significant lower concen- 
trations. We have also studied the infiuence of multi-core 
haloes, i.e., major merger systems which the FOF joins into 
a single system and which other halo-finders would have 
probably singled out each component, in the strength of the 
LFT shape-dynamics correlation and verified that although 
they enhance the correlation, it however persists even ex- 
cluding this subsample of haloes. 

Finally, we find a fraction of the LFT haloes showing 
a somewhat unexpected behaviour, having high sphericities 
and velocity dispersions (comparable or even higher than 
those of more virialized systems), a fact which could be 
attributed to either some late formation-time haloes hav- 
ing transitory spherical shapes and large velocity dispersion, 
during the minimum pericentric passage of a merging pro- 
cess, or/and to a relatively faster virialization process oc- 
curing in some of these haloes. 

The halo morphology- velocity dispersion correlation, al- 
beit weaker, survives the 2D projection, indicating that an 
analogous behaviour should be expected in observational 
group samples. Applying a similar analysis to a "volume- 
limited" subsample of the SDSS DR7 groups of galaxies, 
identified using a friends of friends algorithm, we have also 
found a group shape and dynamics correlation, which is in- 
dependent of group mass (see also Tovmassian & Plionis 
2009). Using a concentration parameter as a proxy of the 
formation time, we indeed find a significant, independent of 
mass, shape-dynamics correlation, in a similar fashion as in 
the dark matter halo case. Less concentrated (late fonuing) 
groups are typically less spherical and have lower velocity 
dispersion than their equal mass more concentrated (early 
forming) counterparts. 
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